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We study theoretically the instability of countersuperflow, i.e., two counterpropagating misci- 
ble superflows, in uniform two-component Bose-Einstein condensates. Countersuperflow instability 
causes mutual friction between the superfluids, causing a momentum exchange between the two 
condensates, when the relative velocity of the counterflow exceeds a critical value. The momentum 
exchange leads to nucleation of vortex rings from characteristic density patterns due to the nonlinear 
development of the instability. Expansion of the vortex rings drastically accelerates the momen- 
tum exchange, leading to a highly nonlinear regime caused by intervortex interaction and vortex 
reconnection between the rings. For a sufficiently large interaction between the two components, 
rapid expansion of the vortex rings causes isotropic turbulence and the global relative motion of the 
two condensates relaxes. The maximum vortex line density in the turbulence is proportional to the 
square of the relative velocity. 

PACS numbers: 03.75.Mn, 67.85. De, 47.27.Cn, 67.25.dk 



I. INTRODUCTION 

Hydrodynamic instability, which causes exotic pat- 
terns or turbulence in the nonlinear development of com- 
plex flow structures, occurs universally throughout na- 
ture from the subatomic to the cosmic scale and has been 
actively studied in many fields such as magnetohydrody- 
namics, plasma physics, elasticity, rheology, and general 
relativity [1]. Such phenomena can occur in quantum 
fluids as well as classical fluids. Quantum turbulence, 
i.e., turbulent states in superfluids, has been thoroughly 
studied mainly in superfluid helium systems and has at- 
tracted considerable attention as an idealized prototype 
of classical turbulence 0, which is the most important 
unsolved problem of classical physics. In quantum flu- 
ids of atomic Bose-Einstein condensates (BECs), there is 
growing interest in the study of the counterparts to hy- 
drodynamic phenomena in classical fluid dynamics, such 
as Kelvin-Helmholtz instability [l, 01, Ray leigh- Taylor 
instability [1, Q , the Benard-von Karman vortex street 
0, and turbulence [8l-[To|. Of greatest importance in 
studying hydrodynamic phenomena in this system is us- 
ing highly developed experimental techniques, that allow 
direct observation of a variety of exotic nonlinear dynam- 
ics caused by macroscopic quantum effects, i.e., superflu- 
idity and quantized vortices. 

Superfluid systems may also exhibit unique instability 
phenomena, which never appear in classical fluid dynam- 
ics. One such phenomenon is the instability of counter- 
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flow, in which two counterpropagating fluid components 
merge into each other. Such a counterflow state is dif- 
flcult to obtain in classical fluid systems because of the 
viscosity between the two components, but is possible in 
superfluid systems where frictionless flow is possible in 
thermal equilibrium. In general, counterflow states be- 
come unstable for a large relative velocity between the 
two components. 

The instability of counterflow is classified into two 
cases. In the first case, the instability occurs between 
a viscid normal fluid component and an inviscid super- 
fluid component. In the second case, the instability is 
between two superfluid components. The former insta- 
bility has in the past been studied in superfluid ^He , 
where the hydrodynamics is usually described using a 
two-fluid model in which the system consists of a vis- 
cous normal fluid component and an inviscid superfluid 
component The counterflow states of the two com- 
ponents are realized by an injected heat current, where 
the two components flow in opposite directions to reduce 
the temperature gradient. When the relative velocity be- 
tween the two components exceeds a critical value, the 
thermal counterflow becomes unstable, developing into 
quantum turbulence with a tangle of quantized vortices. 
In the transition to the turbulent states, remnant vor- 
tices, i.e., quantized vortices attached to the container 
wall, are stretched by the mutual friction between the 
superfluid and normal fluid components, growing into a 
tangle through reconnections with other vortices. Coun- 
terflow turbulence was visualized only recently in exper- 
iments |13l] and is providing an active ground of study in 
quantum turbulence |14| - [l6| . 

It is interesting that counterflow instability causes mu- 
tual friction, i.e., decay of relative motion, even in coun- 
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tersuperflow, which is counterflow between two superflu- 
ids without viscid normal fluids. Recently, we have re- 
ported that the instability of countersuperflow, namely, 
countersupcrflow instability (CSI), can also develop into 
quantum turbulence in gaseous two-component BECs 
[17| . CSI has been studied theoretically in several multi- 
component miscible superfluid systems, including helium 
superfluids [18- 21], mixture BECs of cold atoms [221 . and 
nucleon superfluids in rotating neutron stars [23| . Al- 
though CSI had been studied only theoretically, recently 
Hamner et al. made the observation of CSI experimen- 
tally in gaseous two-component BECs by accelerating the 
two components in opposite directions, utilizing the Zee- 
man shift under a magnetic- field gradient |24| . They 
observed that shocks and dark-bright solitons nucleate 
via the instability in quasi-one-dimensional cigar-shaped 
two-component BECs. Since such a soliton is unstable 
against decay into quantized vortices in two- and three- 
dimensional systems (25j , the CSI can cause vortex nu- 
cleation leading to quantum turbulence even in trapped 
systems with large atomic clouds [iTj . 

This paper presents a detailed analysis of CSI in a 
two-component BEC from linear instability to the non- 
linear development into quantum turbulence in uniform 
systems, focusing on the momentum exchange between 
the two components and the parametric dependence of 
the instability. In Sec. II we review the linear stabil- 
ity of countersuperflow in Bogoliubov analysis. We find 
that the momentum exchange is unique to the counter- 
superflow states. Sections III and IV develop the dis- 
cussion into the nonlinear regime of CSI by numerically 
solving the Gross-Pitaevskii (CP) equations. Section III 
is devoted to density pattern formations and vortex nu- 
cleation in the early stage of the nonlinear dynamics. 
We find that the phase diagram of the linear stability 
characterizes the density patterns depending on the rel- 
ative velocity and the distributions of the resulting vor- 
tex rings typically follow the structures of the patterns. 
Section IV discusses the development into vortex tangles 
from the vortex nucleation. The development of vortex 
lines is similar to that in the thermal counterflow of su- 
perfluid "'He. In our system, the momentum exchange 
between the two components plays the role of mutual 
friction between the normal fluid and superfluid compo- 
nents in thermal counterflow. Even if the total length 
of the initially nucleated vortex line is small, vortex line 
stretching leads to vortex tangles and the momentum 
exchange is completed for a large interaction between 
the two components. The vortex tangles obtained are 
isotropic, in contrast to the vortex tangles under thermal 
counterflow. The results are summarized in Sec. V. 



II. COUNTERSUPERFLOW INSTABILITY 

Let us consider two counterpropagating superfluids 
consisting of two distinguishable particles in an isolated 
uniform system at zero temperature. When the relative 



velocity between the two superfluids exceeds a critical 
value, the countersuperflow becomes unstable. The in- 
stability should then induce a momentum exchange be- 
tween the two superfluids to decay their relative motion. 
Note that CSI is essentially different from Landau insta- 
bility 12], in which the superfluid system is treated in 
the framework of a canonical ensemble. In Landau insta- 
bility, superfluids decay due to friction with the "rigid" 
environment, e.g., the container wall or a normal fluid 
component whose velocity field is frozen to the wall due 
to its viscosity. On the other hand, CSI is purely an 
internal instability with no influence from the external 
environment. Thus, the total energy of the two superflu- 
ids is conserved, keeping its total momentum constant, 
because of Galilean invariance. 

It is physically evident that the momentum exchange 
in CSI cannot initiate a reduction of the relative motion 
of the superfluids as a whole. The momentum exchange 
must arise from a gradual excitation of internal motions, 
i.e., from the appearance of excitations in the superflu- 
ids. In the CSI of two-component BECs, the momentum 
exchange due to the excitations is analytically evaluated 
with the CP and the Bogoliubov-de Gennes (BdG) mod- 
els. In this section, we discuss the linear stability and the 
momentum exchange in the CSI of two-component BECs 
in the BdG models. 



A. Formulations 

We consider a binary mixture of BECs described by the 
condensate wave functions ^j(r,t) = ■^nj(r, t)e*'^J' 
in a mean-fleld approximation at T = K, where the 
index j refers to each component (j — 1,2). The wave 
functions are governed by the coupled CP equations [l^ 
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(1) 



where fc = 1, 2 and nij is the mass of the jth component. 
The coefficient gjk — 2TTlT?ajf^/mjk represents the atomic 
interaction with mjf} = mj^ +"^fc ^ ^^id the s-wave scat- 
tering length Qjk between the jth and kth components. 
Our analysis satisfles the conditions 311322 > 312 and 
gjj > that the two miscible condensates are stable [1^ . 

The wave functions ^'j = in a stationary state are 
written as 



^i{mj'Vj-r — fj,jt)/h 



(2) 



with the velocity Vj and the chemical potential fij of the 
jth component. The countersuperflow is realized with 
Vi ^ V2. We consider a collective excitation above the 
stationary state as ^ j = -I- where we may write 
the excitation wave functions 5'^ ^ with the usual form 



J{mjVyr-lijt)/h 



*}-(3) 
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By linearizing the GP quation [Eq. ([T])] with respect to 
(5^j, we obtain the BdG equations, whose matrix nota- 
tion is 



aMW = eW, 



(4) 



where V is the system volume and we used the particle 
number conservation law of the jth component 



drl*^" + <5vl.,f = V{n, + |u,f + |«,f ) = N,. (10) 



where 



M 



( ht 

giin-i 



giini 



9i2^/nin2 gi2^/nin2 



gi2^Jnin2 gi2^/nTn2 522^^2 

V gi2y/nin2 gi2y/nin2 522"-2 



In an isolated uniform system, the momentum conser- 
vation law requires dJi + 6J2 = 0. Therefore, for the 
momentum exchange, we need the following condition to 
be satisfied: 
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The condition can be evaluated using the quantity 
W'fMW - (W^'MW)* with M = TW^, which is reduced 
to 
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Here Vr — Vi — V2 is the relative velocity, — 
li^q^ /2raj, and e — huj — Vq ■ hq. We may neglect 
the center-of-mass motion of the two condensates Vg = 
{piVi + p2V2)/{pi + P2) with pj = rrijUj, so we set Vj 
and pj to eliminate the velocity of the center of mass. 
The eigenvector W is satisfied with 
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(5) 
(6) 



Because the operator aAi in Eq. (U) is non-Hermitian, 
the eigenvalue e can have an imaginary part. When the 
eigenvalue has an imaginary part, the system is dynami- 
cally unstable. 



B. Momentum exchange 

We investigate how the momentum exchange happens 
via an excitation. We can define the momentum Jj of 
the jth component as 

J, = / dr|(vI/*VvE', - vl^^VvI.*). (7) 

For the steady countersuperflow state, we have 

Jtotal- Jl+J2=A^imiFi+iV2TO2V2, (8) 

with Nj — J dr\'i>j\^. Here we consider that a collective 
excitation of Eq. ([3]) appears in this state. The change 
of momentum of the jth component by the excitation is 
then 



SJ, =V{\u,\^-\v,\^)hq, 



(9) 



{e-e*){\ni\^-\vi\^ + \n2\^-\v2\^) = 0. 



(12) 



, ei) ± ^{el - 4) + lQe\elnin2gl2 with 
= (e^)^ + '^^'jgjj'i^j [2a|- Although e is purely imagi- 



When the eigenvalue e has a complex value e 7^ e*, we 
have [wip — |wip-l-|M2p— |w2p =0 and then the excitation 
is possible satisfying the momentum conservation law. 

Additionally, relative velocity is necessary to exchange 
momentum between the two condensates. The disper- 
sion relation for Vr = is well known and desc ribed 
bye^ = i(e? 

nary when the interaction parameters satisfy the condi- 
tion 311522 < .9i2i thfi excitations then do not cause mo- 
mentum exchange, i.e., SJj = with \uj\ = \vj\, which 
is derived from Eqs. ([5]) and ([6]). Accordingly, the con- 
dition ()lip is satisfied only in dynamic instability with a 
finite relative velocity between the two components. 



C. Dispersion relation and phase diagram 

We can obtain the dispersion relation by solving the 
eigenvalue problem of the BdG Eqs. (HJ, which yields 
a quartic equation in e. Although the eigenvalue prob- 
lem can in principle be solved analytically, the explicit 
expression is generally quite complicated. Here, we con- 
sider the symmetric parameters with tou = 77122 = 
ni = 77.2 = n, and gn — 522 = 5, which are realistic in 
experiments with *^Rb atoms. In this case, we obtain a 
simple form of the dispersion relation 
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- q'%^Vg + q'^l\ 



(13) 



where e' — e/gn 
h/y/mgn, = 

and 7 



q^ with healing length ^ = 
|Vr|/c with sound velocity 



312/5. Here q'^ — + q'1 with 



W ' ^rI^r\ ^'^d l'± — 0- Then, the only important 
parameters in the countersuperflow are = Vr/c and 
7 — 512/5- We have Im e' ^ if the right hand side 
of Eq. (|13p becomes negative. By comparing the square 
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FIG. 1: (Color online) Phase diagram of the countersuperflow 
instability for 7 — 0.9 and 0.1. The gray (purple) region 
shows the unstable region, which shows the inequality (|14|l . 
The vertical axis is the relative velocity — |V)^| — |Vr|/c 
with c = ^/gn/m. The horizontal axis is the wave number 
g|l j_ — g||,xC of the excitations, where we use ^ — h/y/mgn 
and gfi^ + q'l = q'^ with gf, = |q' • Vj^/l^Al and > 0. The 
color bars show the magnitude of the imaginary part of the 
eigenvalues e' — e/gn. The dashed lines represents = 4.71 
[(a) and (b)] and 2.36 [(c) and (d)] and the solid lines show 
the critical velocities V± — 2-^1 ± 7. The lower four figures 
show the cross-section surfaces of = 2.36 and 4.71 in the 
phase diagrams. 



of the fourth term to that of the sum of the other three 
terms on the right-hand side of Eq. (fT3|) . we obtain the 
condition of the CSI 



Figure [T] shows this condition with 7 = 0.9 and 0.1, 
where the gray (purple) zone shows that the countersu- 
perflow is unstable. When is small, no mode is un- 
stable. However, some modes become unstable when 
exceeds the critical value V!_ with 



2v/(T±^, 



which is based on Eq. (|14l) . 

For small relative velocity with V!_ < V^^ < V^, the 
unstable region is broadly distributed to the lower-wave- 
number region [Fig. [I](c)]. When increases to suffi- 
ciently large relative velocity with > V_l, the unsta- 
ble modes with larger values of |Im e'\ are distributed to 
the higher- wave-number region in the wave-number space 
((7||,q^) and the cross section forms crescents [Figs. [Tfa) 
andlllb)]. 



III. VORTEX NUCLEATION 

In Sec. II we demonstrated the linear stability analy- 
sis of counterpropagating two-component BECs. In the 
linear stage of the instability, the unstable modes are 
amplified exponentially with time and momentum is ex- 
changed between the two components. It is expected that 
the wave numbers of the unstable modes with larger val- 
ues of |Im e'l could affect the density patterns emerging 
after the amplification in the nonlinear regime of the in- 
stability. Hamner et al. have observed experimentally 
that the CSI causes characteristic solitary waves in the 
nonlinear stage in a quasi-one-dimensional system of two- 
component BECs [24| . In a three-dimensional system, 
the density patterns lead to vortex nucleation in the early 
stage of the nonlinear development [r?j . 

In this section we discuss the nonlinear dynamics of 
CSI until vortex nucleation in a three-dimensional sys- 
tem by numerically solving the coupled GP equations 
[Eq. ([T|)]. The numerical simulations were done in a 
three-dimensional system, which was subject to a peri- 
odic boundary condition. The initial state is the sta- 
tionary state of Eg. (|p with small white noise to trig- 
ger the instability (27[, and we set the parameters as 
TOi = m2 = m, 511 = 522 = g, n-i = n2 = n, V{ = -V2, 
and jij = myJ/8 + {g + gi2)n. Because the parameters of 
two components are symmetric, both components follow 
similar scenarios in the nonlinear dynamics. 

The vortex nucleation is generally classified into two 
cases, > and < < V^, according to 
the configuration of the region of the unstable modes 
in the phase diagram (Fig. [1]). First, we demonstrate 
the vortex nucleation dynamics for a large relative veloc- 
ity > V+- Figures [2I^a)-[2{d) show the typical vortex 
nucleation dynamics with — 4.71 > along the x 
axis, l^'jl^/n = 1, and 7 — 0.9. In the early stage of 
the nonlinear development after the exponential ampli- 
fication of the unstable modes, disk-shaped low-density 
regions appear in both components, which face in the di- 
rection parallel to the relative velocity Vr [Fig. [IJb)]. A 
tiny vortex ring is nucleated in the disk region [Fig. [Sfc)] , 
which immediately grows to a vortex ring with a certain 
size [Fig. (Hd)]. Since the vortex nucleation results from 
the amplification of the unstable modes, decreasing the 
relative motion, the nucleated vortex rings in the jth 
component must carry the momentum in the opposite 



(15) direction to the initial velocity Vj, where the direction of 
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Component 1 (vortex core and density iso.'^urface) 




Compononi I and 2 (voncx core only) 




(a) t' = (b) t' - 12.2 Cc) i' = 12.8 (d) t' = 13.,^ (e) t' = 13.8 (0 = 26.0 

FIG. 2: (Color online) Nonlinear dynamics of vortex cores in countersuperflow instability. Curved lines in the top and middle 
panels show the vortex cores of the first and second components, respectively. The vortex cores of both component are plotted 
together in the bottom panels. The curved surfaces in top top and middle panels represent the density isosurfaces of the first 
component with I^Pil^/n = 0.1 and the second component with |*I'2|^/n — 0.1, respectively. First, the instability causes disk- 
shaped low-density regions that face in the direction parallel to the initial relative velocity Vr \\ x with a unit vector x along 
the X axis (b). Then small vortex rings appear inside the low-density disks (c) and the disk then immediately transforms into a 
torus-shape like a usual vortex ring (d) . The nucleated vortex rings gradually expand due to momentum exchange between the 
two components. The vortices then start to reconnect with each other (e). Frequent reconnections then cause vortex tangles 
(f). The time t is normalized as t = tr with r — h/gn. The box size is V = (8^)"^ with ^ — h/^/mgn. 



the phase gradient Vc^t inside the vortex ring is opposite 
that outside the ring 

Figure [3] shows a detailed typical dynamics of the vor- 
tex ring nucleatioin in the first component. First, a disk- 
shaped low-density region appears and the phase gradient 
inside the disk region becomes steep [Fig. [31[a)]. Then a 
vortex ring appears from the zero-density point inside the 
disk region, where the radius of the vortex ring is smaller 
than that of the disk region [Fig. EJb)]. The radius of the 
ring grows immediately into that of the disk region and 
the geometry of the region transforms from a disk into a 
torus, like a usual vortex ring [Fig. El^c)]. Therefore, the 
size of the resulting vortex ring is roughly characterized 
by that of the disk-shaped low-density region. 

Under this consideration, we may estimate the dis- 
tribution of vortex rings from the structure of the den- 
sity patterns emerging in the early stage. The structure 



may be determined by the wave numbers of the unstable 
modes with large values of jlm s'\ in the phase diagram 
(Fig. [T|). In the limit of large relative velocity V^^^ V^, 
the unstable region described with the inequality in Eq. 
(fT4| is reduced to the circle 

iiii-lykr+<ii = lv^'- (16) 

As a result, the wave numbers of the unstable modes 
are characterized by q'^^ < and < for suffi- 

ciently large relative velocity. Figure |4l[a) shows a typical 
density pattern in a cross-section plane along the relative 
velocity (top panels) and its vortex plot (bottom panels) . 
Most of the low-density disks nucleate vortex rings. Thus 
the periodicities ~ Vr/k of the density pattern charac- 
terize the vortex ring size and the interval between the 
vortex rings, where we used n = 2-Kh/m. Therefore, for 
Vjzj > V!^, the vortex line density just after the nucleation 



event may be roughly estimated to be ~ V^/k^, which is 
independent of 7. 




FIG. 3: (Color online) Vortex ring nucleation from a low- 
density disk in countersuperflow instability. Top panels show 
the time development of the vortex core and the density iso- 
surface of the first component with |^i|^/n = O.f. Middle 
panels represents plots of the density l^'ij^/n and phase phii 
in the cross section cutting the ring in half. Shown in the bot- 
tom panels is a schematic diagram of the vortex nucleation. 
The shaded areas and solid curves show the low-density re- 
gions and the 7r/2-spaced isolines of phase, respectively. The 
density patterns due to the CSI induce the disk-shaped low- 
density region. The phase gradient is steep inside the low- 
density disk. The disk decays into a vortex-antivortex pair 
(a vortex ring) whose size corresponds to the disk size. After 
the vortex nucleation, the phase gradient becomes smaller and 
the superflow is reduced between the pair (inside the ring). 

However, the situation is different when 7 or Vr be- 
comes small. When the interaction between the two com- 
ponents is weak for small 7, the density pattern is highly 
disturbed by the time the vortex rings are nucleated [Fig. 
IHb) (top)]. As a result, the number of nucleated rings in 
Fig. m^b) is smaller than that in Fig. SJa), although the 
rings have the same characteristic size ~ Vr/k. 

If the relative velocity is small, e.g., < < V^, 
since the unstable modes with large |Im e'| are dis- 
tributed to low values of q'j^ [Fig. [He)], the periodicity 
of the density pattern perpendicular to the relative ve- 
locity becomes larger than that parallel to the velocity 
[Fig. mjc) (top)]. Thus the radius of nucleated vortex 
rings is larger than the interval between the rings along 
the relative velocity [Fig. |3Jc) (bottom)], where the vor- 
tex line density becomes small compared to V^/k'^. If 
the perpendicular periodicity of the density pattern be- 
comes comparable to the system size, the CSI does not 
make vortex rings, but rather makes vortex lines across 
the system or distorted stripe patterns, as observed by 



Hamner et al. j24| . 
(a) Y=0.9, Vr'=4.71 (b) 7=0.1 , Vr'=4,7 1 (c) 7=0.9, Vr'=2.36 
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FIG. 4: (Color online) (Top panels) Density pattern of the 
first component in a cross section along the initial relative ve- 
locity. (Bottom panels) Plots of cores of the nucleated vortices 
of the first component, (a) and (b) For large relative velocity 
Vfi > V4. , the low-density region is small in the direction per- 
pendicular to the relative velocity and small vortex rings are 
nucleated, (c) For small relative velocity V4. > Vli > VL, the 
density pattern is stretched in the direction perpendicular to 
the relative velocity and a small number of large vortex rings 
are nucleated. The box size is V = (f6^)^. 



IV. TRANSITION TO VORTEX TANGLE 

Since the momentum transported by a vortex ring in- 
creases as its radius increases ^llj . further momentum 
exchange expands the nucleated vortex rings without 
changing their direction. As the vortex rings expand, 
the intervortex interaction between neighboring vortex 
rings becomes important, leading to highly nonlinear dy- 
namics. 

The typical dynamics after vortex nucleation are 
shown in Figs. [2Ud)-[2jf). At the beginning of the ex- 
pansion, the size of the nucleated vortex rings gradu- 
ally becomes large without changing their direction [Fig. 
I^d)]. As the vortex rings expand further, they are more 
distorted by the interaction and reconnection with neigh- 
boring rings [Fig. [D^e)]. The vortex distortions are more 
affected by the vortices in the same component than by 
those in the different components. This is because the 
force between vortices in the same component is stronger 
than that in the different components for large distances 
[28| . In addition, vortices in the same component can 
make the usual reconnections [llj , while those in the dif- 
ferent components just cross each other even if they col- 
lide. The distortion leads to a tangle of quantized vortices 
and then both components follow the same sequence to 
form binary quantum turbulence [Fig. [^f )] 
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This scenario is similar to that of the transition to 
quantum turbulence in thermal counterflow 11|. In 



(a) 7=0.9. Vr'=4.71 (b) y=0.1.Vh'=4.71 (c) y=0.9,VR'=2.36 



the transition dynamics to thermal counterflow turbu- 
lence, the vortex lines attached to the container wall are 
dragged and stretched by mutual friction, which causes a 
momentum exchange between the superfluid and normal 
fluid components. In our system of CSI, the momentum 
exchange between the two components plays the role of 
mutual friction in the thermal counterflow. The mutual 
friction is microscopically understood as the reaction of 
scattering of the excitations, of which the normal fluid 
component consists, at the core of the quantized vor- 
tices. However, this picture cannot be simply applied to 
the "mutual friction" in two-component BECs because 
neither of the two components is a normal component. 
A microscopic understanding of the "mutual friction" is 
an open problem for future work. 

There are two evident differences between CSI and 
thermal counterflow instability. One difference is the 
vortex nucleation mechanism. In CSI, the vortices are 
intrinsically nucleated in the bulk, as discussed above, 
whereas in thermal counterflow the vortices are prepared 
on the boundary of the system before the instability. 
The second difference involves the relative motion be- 
tween the two components. Thermal counterflow turbu- 
lence has been traditionally studied in states in which 
the superfluid component is turbulent but the normal 
fluid component is laminar sustained by an externally 
driving temperature gradient. For CSI, in contrast, the 
relative motion decays through the momentum exchange 
between the two components and both components be- 
come turbulent, keeping its total momentum conserved 
in an isolated uniform system [29j . 

In order to obtain further insight into CSI, we inves- 
tigate how the relative motion decays, stretching vor- 
tices throughout the nonlinear dynamics. The typi- 
cal time development of the mean superfluid velocity 
Vx,j = Jx,j /iTT'jNj with the momentum J^.j of the jth 
component in the x direction is plotted together with 
the total vortex line density I in Fig. ^a) (middle) . The 
rate of the momentum exchange drastically increases af- 
ter the vortex nucleation and decreases when vortex re- 
connection takes place frequently. This implies that the 
expansion of vortex rings is the most efficient mechanism 
for momentum exchange in CSI. After vortex tangles are 
realized via frequent vortex reconnection, the momentum 
exchange moderately continues even when the averaged 
relative velocity becomes lower than the critical relative 
velocity, i.e., \vx^i — Vx,2\ < and then the momentum 
exchange is almost completed and the relative motion 
vanishes globally. 

In the momentum exchange process, the kinetic en- 
ergy due to the global relative motion with v^^i ^ 
is mainly consumed in driving the internal motion, e.g., 
nucleation and stretching of vortices in each component. 
This is why the growth rate of the vortex line density fol- 
lows the rate of momentum exchange in Fig. [SJa) (mid- 
dle), e.g., the growth rate of the vortex line decreases 










0.3 




2 




11.15 (1 




-2 










0.4 



0.2 



30 n 100 200 



0.6 

0,5 

Oil 





0.6 




0,5 




0.3 




1) 



20 



30 



100 



150 



200 




FIG. 5: (Color online) (Top panels) Plots of vortex cores of 
the first component at the time when the vortex line density 
is maximum. The box size is V = (160^ (Middle panels) 
Time evolution of the total vortex line density I' — of 
both components and the normalized mean superfluid veloc- 
ity v'^ j^^/t = Jx.j /rrijNj with the momentum Jx,j of the jth 
component in the x direction. (Bottom panels) Time evo- 
lution of the parameter Q\\/Q±. The momentum exchange 
happens drastically in the expansion of the nucleated vortex 
rings. The rate of the momentum exchange is reduced when 
vortex reconnection occurs frequently in the vortex tangles. 
Eventually, the momentum exchange is completed, leading to 
an isotropic tangle. For small 7 = 0.1 (b), the vortex rings 
expand less and thus the momentum exchange is completed 
without an isotropic tangle. Even if the vortex line density 
just after the vortex nucleation event is small for small rela- 
tive velocity V4. > > ^- with large 7 = 0.9, the vortices 
are well stretched and repeat reconnections to develop into a 
dilute isotropic tangle, where the momentum exchange is also 
completed. 



when the rate of momentum exchange becomes small af- 
ter frequent vortex reconnection. Generally speaking, the 
vortex line length is reduced by emitting phonons via vor- 
tex reconnection pd| and hence the vortex line density 
eventually stops growing when the reduction due to vor- 
tex reconnection exceeds the growth effect by momentum 
exchange. 

The maximum value Zmax of the vortex line density 
should increase with the initial relative velocity Vr since 
the initial kinetic energy is consumed in stretching vor- 
tex lines. Figure [5] shows the numerical result of the 
maximum vortex line density ^max dependence on Vr 
for 7 — 0.9. The vortex line density ^max is propor- 
tional to V^. This behavior of ^max is phenomenologi- 
cally explained from the energetics considering the ten- 
sion of a vortex line. To apply this consideration to the 
CSI turbulence, we assume that all of the kinetic en- 
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max 




FIG. 6: (Color online) Maximum total vortex line density 
^max = ^max^^ in the CSI turbulence with 7 = 0.9 as a function 
of , which is obtained by making an ensemble average of 
five numerical results with different initial noise. 



ergy due to the global relative motion is used for stretch- 
ing vortices. Then the kinetic-energy density in the 
initial state is equal to aZmax, where the proportional- 
ity factor a represents the effective tension coefficient 
of the vortex line. Since the energy density vq is writ- 



ten as Vf) 



P1+P2 



= mn/2, we obtain 



with Vg = and 



(17) 



where we have used a oc nh^/m as an analog to the 
tension of an isolated vortex line. The numerical result 
of Fig. inimay imply that a is not sensitive to I, although, 
in general, a depends logarithmically on the value a/^„ 
with the mean intervortex distance a = l''^/'^ and the 
vortex core radius ^t, [Tlj . 

The usual thermal counterflow turbulence is 
anisotropic in steady states because relative motion 
remains between the superfluid and normal fluid com- 
ponents, sustained by an externally driving temperature 
gradient It is expected that CSI turbulence can be 
isotropic since the relative motion vanishes in the final 
stage. We investigate the anisotropy of CSI turbulence in 
a manner similar to the analysis of thermal counterflow 
turbulence by introducing the anisotropic parameters 

m 



Q 



1 



1 



where we have used the system volume V and the effec- 
tive vorticity u) = (w^, Wj^, tiJ^)"^ = V x v with the total 



mass-current velocity D(r) = {piVi + P2V2) / {pi+ P2) and 
Vj{r) = ti/mj'V(j)j{r). Analogous to thermal counter- 
flow turbulence, we have Q\\/Q± = 1/2 for isotropic CSI 
turbulence. Conversely, if the tangle consists of curves 
lying only in planes perpendicular to Vr, then we have 

Q\\/Q±^o. 

We show the typical time development of Q\\/Q± from 
the vortex nucleation in Fig. [5] (a) (bottom). Q\\/Q± is 
almost zero just after the vortex nucleation because the 
vorticity lo of the nucleated vortex ring is perpendicular 
to the relative velocity. After the vortex rings start to 
interact with other vortices and get distorted, the vor- 
tex tangle grows to be isotropic Q\\/Q± ~ 1/2 as the 
momentum exchange is completed Vx,i = Vx.2 = 0. 

The above analysis is not applicable to the smaller 7 
case, where the nucleated vortex rings expand little and 
reconnect less with the other vortex rings [Fig. [5] (b)]. 
In this case, vortex lines tend to lie perpendicular to the 
relative velocity with Q\\/Qi_ <^ 1/2, where there remain 
vortex rings facing parallel to the relative velocity, as 
seen in Fig. [5] (b) (top) . Then the momentum exchange 
is not completed and slows even when the averaged rela- 
tive velocity is still much larger than the critical relative 
velocity: {v^^i - Vx,2\ > 

In contrast, the relative velocity has less influence on 
the nonlinear dynamics. For small Vn with 7 = 0.9, the 
vortex rings expand readily and repeat reconnections to 
develop into isotropic turbulence without global relative 
motion of the two components [Fig. [5] (c)]. It is inter- 
esting that the vortex tangles are anisotropic in Fig. [5] 
(b) but isotropic in Fig. [5](c), although the values of the 
maximum vortex line density are comparable for both 
cases. This shows that large 7 (< 1) rather than large 
Vji is crucial for realizing an isotropic tangle. 



V. SUMMARY 

We studied theoretically the instability of countcrsu- 
perflow in uniform miscible two-component BECs. It was 
shown that the dynamic instability induces momentum 
exchange between the two components with a finite rel- 
ative velocity in the Bogoliubov analysis. The nonlinear 
development of the instability was analyzed by numeri- 
cally solving the coupled CP equations. The distribution 
of the nucleated vortex rings was qualitatively evaluated 
from the instability phase diagram. The size of the nucle- 
ated vortex rings increases as the initial relative velocity 
Vr decreases. The nonlinear dynamics after vortex nu- 
cleation are similar to those of the turbulence transition 
in thermal counterfiow of superfiuid "^He. In contrast to 
the anisotropy of thermal counterfiow turbulence, CSI 
causes isotropic turbulence by completing the momen- 
tum exchange. The maximum vortex line density of the 
turbulence is proportional to (Vr/k)^. Even if the total 
line length of the nucleated vortices is small for small Vr, 
isotropic turbulence is realized by vortex stretching due 
to mutual friction for a sufficiently large intercomponent 
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interaction. 

In this work we considered uniform density conden- 
sates. However, density is inhomogeneous in the experi- 
mental reahzations of BEC in dilute atomic clouds. Our 
results may still be applied to the case in which the size of 
an atomic cloud is large enough compared to the density 
pattern periodicity ~ Vr/k parallel and perpendicular 
to the relative velocity Vr with large Vr > V+. Since 
the critical relative velocity decreases with the conden- 
sate density, the critical relative velocity at the surface of 
the condensates is smaller than that at the center in the 
local density approximation. Thus the instability can oc- 
cur from the surface before the relative velocity exceeds 
the critical relative velocity at the center. Then vortices 
would be nucleated not in the bulk but from the surface 
after the amplification of the unstable surface modes. 
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